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A Cyclic Model of the Universe 


Paul J. Steinhardt’* and Neil Turok? 


We propose a cosmological model in which the universe undergoes an endless 
sequence of cosmic epochs that begin with a “bang” and end in a “crunch.” 
Temperature and density at the transition remain finite. Instead of having an 
inflationary epoch, each cycle includes a period of slow accelerated expansion 
(as recently observed) followed by contraction that produces the homogeneity, 
flatness, and energy needed to begin the next cycle. 


The current standard model of cosmology 
(1-4) combines the original big bang mod- 
el and the inflationary scenario. Inflation, a 
brief period (10° °° s) of very rapid cosmic 
acceleration shortly after the big bang, can 
explain the homogeneity and isotropy of 
the universe on large scales (>100 Mpc), 
its spatial flatness, the distribution of gal- 
axies, and the spatial fluctuations in the 
cosmic microwave background. However, 
the standard model has some cracks. The 
recent discoveries of cosmic acceleration 
and gravitationally self-repulsive dark en- 
ergy (5—8) were not predicted and have no 
particular role in the standard model (1—3). 
Furthermore, the standard model does not 
explain the “beginning of time,” the initial 
conditions of the universe, or what will 
happen in the long-term future. 

Here, we present a cosmological model 
with an endless sequence of cycles of ex- 
pansion and contraction. By definition, 
there is neither a beginning nor end of time, 
nor is there a need to define initial condi- 
tions. In addition, we explain the role of 
dark matter and generate the homogeneity, 
flatness, and density fluctuations without 
invoking inflation. 

The cyclic aspect of the model is reminis- 
cent of oscillatory models introduced in the 
1930s based on a closed universe that under- 
goes a sequence of expansions, contractions, 
and bounces. The oscillatory models were 
constrained by having to pass through a sin- 
gularity in which the energy and temperature 
diverge. Furthermore, as pointed out by Tol- 
man (9, 10), entropy produced during one 
cycle would add to the entropy produced in 
the next, causing each cycle to be longer than 
the one before it. Extrapolating backward in 
time, the universe would have to have origi- 
nated at some finite time in the past so that 
the problem of explaining the “beginning of 
time” remains. Furthermore, today we know 
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that the universe is flat, rather than closed, 
based on measurements of the cosmic micro- 
wave background anisotropy and large-scale 
structure. 

In our cyclic model, the universe is infi- 
nite and flat, rather than finite and closed as 
in the oscillatory models. We introduce a 
negative potential energy rather than spatial 
curvature to cause the reversal from expan- 
sion to contraction. Before the reversal, 
though, the universe undergoes the usual pe- 
riod of radiation and matter domination, fol- 
lowed by a long period of accelerated expan- 
sion [presumably the acceleration that has 
been recently detected (5—8)]. The accelerat- 
ed expansion, caused by dark energy, is nec- 
essary to dilute the entropy, black holes, and 
other debris produced in the previous cycle so 
that the universe is returned to its original 
pristine vacuum state before it begins to con- 
tract, bounce, and begin a cycle anew. 


Essential Ingredients 


As in inflationary cosmology, the cyclic sce- 
nario can be described in terms of the evolu- 
tion of a scalar field ¢ along a potential V( ) 
in a four-dimensional (4D) quantum field 
theory. The essential differences are the form 
of the potential and the couplings between the 
scalar field, matter, and radiation. 

The analysis of the cyclic model follows 
from the action S that describes gravity, the 
scalar field @, and the matter and radiation 
fluids: 
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where g is the determinant of the metric 
Zp» G is Newton’s constant, and R is the 
Ricci scalar. The coupling B( ) between b 
and the matter (p,,) and radiation (ppg) den- 
sities is crucial because it causes the den- 
sities to remain finite at the big crunch/big 
bang transition. 

The line element for a flat, homogeneous 
universe is —dit? + a?dx*, where a is the 
Robertson-Walker scale factor. The equa- 
tions of motion following from Eq. 1 are 


where a dot denotes a derivative with respect 
to t and H = /a is the Hubble parameter. The 
equation of motion for ¢ is 
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and the fluid equation of motion for matter 
(M) or radiation (R) is 
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where â = aßB(ẹ), p; is the pressure of the 
fluid component with energy density p, and 
“b” is d/dẹọ. The implicit assumption is that 
matter and radiation couple to B()g,,, 
(with scale factor â) rather than the Einstein 
metric g alone (or the scale factor a). Note 
that the radiation term in Eq. 1 is actually 
independent of p (because pp * G *) so only 
Py, enters the p equation of motion. 

We assume the potential V(b) has the 
following features, as illustrated in Fig. 1: (i) 
V(b ) must approach zero rapidly as þ—> —%; 
(ii) the potential must be negative for inter- 
mediate ; and (iii) as increases, the po- 
tential must rise to a shallow plateau with a 
positive value V,. An example of a potential 
with these properties is 


Vib) = Vo -e**) Fb) (6) 


where from this point onward we adopt units 
in which 87G = 1. F(d) is a function we 
introduce to ensure that Vib)— 0 as b> 
—, Without loss of generality, we take F( ) 
to be nearly unity for ¢ to the right of poten- 
tial minimum. The detailed manner in which 
it tends to zero at smaller ¢ is not crucial for 
the main predictions of the cyclic model. A 
quantitative analysis of this model potential 
(11) shows that a realistic cosmology can be 
obtained by choosing c = 10 and V, equal to 
today’s dark-energy density (~6 X 10730 
g/cm?) in Eq. 6. 

We have already mentioned that the cou- 
pling B(¢) is chosen so that â and, thus, the 
matter and radiation density are finite at a = 
0. For example, we consider B() ~ e~?” V° 
as þ—> —%, The presence of B(d) and the 
consequent coupling of to nonrelativistic 
matter represent a modification of Einstein’s 
theory of general relativity. Because the sca- 
lar field ọ evolves by an exponentially small 
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amount between nucleosynthesis (tf ~ 1 s) 
and today (t ~ 10!7 s), the deviations from 
standard general relativity are small enough 
to easily satisfy all current cosmological con- 
straints (77). However, the coupling of matter 
to b produces other potentially measurable 
effects, including a fifth force that violates 
the equivalence principle. Provided (In B) 
< 1073, for today’s value of ¢, these vio- 
lations are too small to be detected (71-13). 
We shall assume this to be the case. Hence, 
the deviations from general relativity are neg- 
ligible today. 

The final crucial ingredient in the cyclic 
model is a matching rule that determines how 
to pass from the big crunch to the big bang. 
The transition occurs as þ—> —© and then 
rebounds toward positive . Motivated again 
by string theory (see below), we propose that 
some small fraction of the -field kinetic 
energy is converted to matter and radiation. 
The matching rule amounts to 


be > -0 + bee 7) 


where x is a parameter measuring the effi- 
ciency of production of radiation at the 
bounce. Both sides of this relation are finite at 
the bounce. 


Stringy Motivation 


From the perspective of 4D quantum field 
theory, the introduction of a scalar field, a 
potential, and the couplings to matter in the 
cyclic model is no more arbitrary or tuned 
than the requirements for the inflationary 
models. However, the ingredients for the cy- 
clic model are also strongly motivated by 
string theory and M-theory. This connection 
ties our scenario into the leading approach to 
fundamental physics and quantum gravity. 
The connection should not be overempha- 
sized. String theory is not yet proven and 
quantum gravity effects are unimportant for 
describing cosmology at wavelengths much 
longer than a Planck length (10-3? cm). If the 
reader prefers, the connection to string theory 
can be ignored. On the other hand, we find 
the connection useful because it provides a 
natural geometric interpretation for the sce- 
nario. Hence, we briefly describe the relation. 

According to M-theory, the universe con- 
sists of a 4D “bulk” space bounded by two 
3D domain walls, known as “branes” (short 
for membranes), one with positive and the 
other with negative tension (74-17). The 
branes are free to move along the extra spatial 
dimension, so that they may approach and 
collide. The fundamental theory is formulat- 
ed in 10 spatial dimensions, but 6 of the 
dimensions are compactified on a Calabi-Yau 
manifold, which for our purposes can be 
treated as fixed, and therefore ignored. Grav- 
ity acts throughout the 5D spacetime, but 
particles of our visible universe are con- 
strained to move along one of the branes, 
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sometimes called the visible brane. Particles 
on the other brane interact only through grav- 
ity with matter on the visible brane and hence 
behave like dark matter. 

The scalar field b is naturally identified 
with the field that determines the distance 
between branes. The potential V() is the 
interbrane potential caused by nonperturba- 
tive virtual exchange of membranes between 
the boundaries. The interbrane force is what 
causes the branes to repeatedly collide and 
bounce. At large separation (corresponding to 
large o ), the force between the branes should 
become small, consistent with the flat plateau 
shown in Fig. 1. Collision corresponds to b= 
—, But the string coupling g, x e%, with 
y > 0, so g, approaches zero in this limit 
(18). Nonperturbative effects vanish faster 
than any power of g, for example, as or 
e "85, or e7 "8s, accounting for the prefactor 
F(d) in Eq. 6. 

The coupling B(¢) also has a natural in- 
terpretation in the brane picture. Particles 
reside on the branes, which are embedded in 
an extra dimension whose size and warp are 
determined by B. The effective scale factor 
on the branes is â = aB(¢), not a, and å is 
finite at the big crunch or big bang. The 
function B(¢) is in general different for the 
two branes (due to the warp factor) and for 
different reductions of M-theory. However, 
the standard Kaluza-Klein behavior B(b) ~ 
e */V® as b—> — œ is universal, because the 
warp factor becomes irrelevant as the branes 
approach one another (//, 18). 

Most importantly, the brane-world pro- 
vides a natural resolution of the cosmic sin- 
gularity (71, 18). One might say that the big 
crunch is an illusion, because the scale factors 
on the branes (— 4) are perfectly finite there. 
That is why the matter and radiation densi- 
ties, and the Riemannian curvature on the 
branes, are finite. The only respect in which 
the big crunch is singular is that the one extra 
dimension separating the two branes momen- 
tarily disappears. Our scenario is built on the 


Fig 1. Schematic plot of the po- 
tential V(þ) as a function of the 
field p. In M theory, ¢ determines 
the distance between branes, 
where þ—> —~ as the branes col- 
lide. We define b to be zero where 
V(¢ ) crosses zero and, therefore, 
is positive when the branes are at 
their maximal separation. Far to 
the right, the potential asymp- 
totes to Po the current value of 
the quintessence (dark energy) 
density. The solid circles represent 
the dark energy—dominated stage; 
the gray circles represent the con- 
tracting phase during which densi- 
ty fluctuations are generated; the 
open circles represent the phase 


hypothesis (/9) that the branes separate after 
collision, so the extra dimension immediately 
reappears. This process cannot be completely 
smooth, because the disappearance of the 
extra dimension is nonadiabatic and leads to 
particle production; that is, the brane colli- 
sion is partially inelastic. Preliminary calcu- 
lations of this effect are encouraging, because 
they indicate that a finite density of particles 
is produced (20). The matching condition, 
Eq. 7, parameterizes this effect. Ultimately, a 
well-controlled string-theoretic calculation 
(11, 18, 20) should determine the value of x. 


Dark Energy and the Cyclic Model 


The role of dark energy in the cyclic scenario 
is novel. In the standard big bang and infla- 
tionary models, the recently discovered dark 
energy and cosmic acceleration (5—8) are an 
unexpected surprise with no obvious expla- 
nation. In the cyclic scenario, however, not 
only is the source of dark energy explained, 
but the dark energy and its associated cosmic 
acceleration are actually crucial to the con- 
sistency of the model; namely, the associated 
exponential expansion suppresses density 
perturbations and dilutes entropy, matter, and 
black holes to negligible levels. By periodi- 
cally restoring the universe to an empty, 
smooth state, the acceleration causes the cy- 
clic solution to be a stable attractor. 

Right after a big bang, the scalar field ¢ is 
increasing rapidly. However, its motion is 
damped by the expansion of the universe, and 
& essentially comes to rest in the radiation- 
dominated phase (stage | in Fig. 1). There- 
after it remains nearly fixed until the dark 
energy begins to dominate and cosmic accel- 
eration commences. The positive potential 
energy density at the current value of & acts 
as a form of quintessence (2/, 22), a time- 
varying energy component with negative 
pressure that causes the present-day acceler- 
ated expansion. This choice entails tuning Vo, 
but it is the same degree of tuning required in 
any cosmological model (including inflation) 


when the scalar kinetic energy dominates; and the broken circle represents the stage when the 
universe is radiation dominated. Further details of the sequence of stages are described in the text. 
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to explain the recent observations of cosmic 
acceleration (5-8). In this case, because the 
dark energy serves several purposes, the sin- 
gle tuning resolves several problems at once. 
The cosmic acceleration is nearly 100 or- 
ders of magnitude smaller than considered in 
inflationary cosmology. Nevertheless, if sus- 
tained for hundreds of e-folds (trillions of 
years) or more, the cosmic acceleration can 
flatten the universe and dilute the entropy, 
black holes, and other debris (neutron stars, 
neutrinos, and so forth) created over the pre- 
ceding cycle, overcoming the obstacle that 
has blocked previous attempts at a cyclic 
universe. In this picture, we are presently 
about 14 billion years into the current cycle, 
and have just begun the trillions years of 
cosmic acceleration. After this amount of 
accelerated expansion, the number of parti- 
cles in the universe may be suppressed to less 
than one per Hubble volume before the cos- 
mic acceleration ends. Ultimately, the scalar 
field begins to roll back toward —%, driving 
the potential to zero. The scalar field ọ is thus 
the source of the currently observed acceler- 
ation, the reason why the universe is homo- 
geneous, isotropic, and flat before the big 
crunch, and the root cause for the universe 
reversing from expansion to contraction. 


A Brief Tour of the Cyclic Universe 


Putting together the various concepts that have 
been introduced, we can now present the se- 
quence of events in each cycle beginning from 
the present epoch, stage | in Fig. 1. The uni- 
verse has completed radiation- and matter-dom- 
inated epochs during which ọ is nearly fixed. 
We are presently at the time when its potential 
energy begins to dominate, ushering in a period 
of slow cosmic acceleration lasting trillions of 
years or more, in which the matter, radiation, 
and black holes are diluted away and a smooth, 
empty, flat universe results. Very slowly the 
slope in the potential causes ọ to roll in the 
negative direction, as indicated in stage 2. Cos- 
mic acceleration continues until the field nears 
the point of zero potential energy, stage 3. The 
universe is dominated by the kinetic energy of 
Q, but expansion causes this to be damped. 
Eventually, the total energy (kinetic plus nega- 
tive potential) reaches zero. From Eq. 2, the 
Hubble parameter is zero and the universe is 
momentarily static. From Eq. 3, ä < 0, so that 
a begins to contract. So long as a is nearly 
static, the universe satisfies the ekpyrotic con- 
ditions for creating a scale-invariant spectrum 
of density perturbations (79, 23). As the field 
continues to roll toward —%, a contracts and the 
kinetic energy of the scalar field grows; that is, 
gravitational energy is converted to scalar field 
kinetic energy during this part of the cycle. 
Hence, the field races past the minimum of the 
potential and off to —%, with kinetic energy 
becoming increasingly dominant as the bounce 
nears (stage 5). The scalar field diverges as a 
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tends to zero. After the bounce, radiation is 
generated and the universe is expanding. At 
first, scalar kinetic energy density (x 1/a®) 
dominates over the radiation (« 1/a*) (stage 6). 
Soon after, however, the universe becomes ra- 
diation dominated (stage 7). The motion of ¢ is 
rapidly damped away, so that it remains close to 
its maximal value for the rest of the standard 
big bang evolution (the next 15 billion years). 
Then, the scalar field potential energy begins to 
dominate, and the field rolls towards —%, 
where the next big crunch occurs and the cycle 
begins anew. 


Obtaining Scale-Invariant 
Perturbations 


One of the most compelling successes of infla- 
tionary theory was to obtain a nearly scale- 
invariant spectrum of density fluctuations that 
can seed large-scale structure (4, 24-26). Here, 
the same feat is achieved using different phys- 
ics during an ultraslow contraction phase (stage 
3 in Fig. 1) (19, 23). In inflation, the density 
fluctuations are created by very rapid expan- 
sion, causing fluctuations on microscopic scales 
to be stretched to macroscopic scales (4, 24— 
26). In the cyclic model, the fluctuations are 
generated during a quasistatic, contracting uni- 
verse where gravity plays no significant role 
(19). Simply because the potential V(d ) is de- 
creasing more and more rapidly, quantum fluc- 
tuations in ¢ are amplified as the field evolves 
downhill (79, 27, 28). Instabilities in long- 
wavelength modes occur sooner than those in 
short-wavelength modes, thereby amplifying 
long-wavelength power and, curiously, nearly 
exactly mimicking the inflationary effect. The 
nearly scale-invariant spectrum of fluctuations 
in ọ created during the contracting phase trans- 
form into a nearly scale-invariant spectrum of 
density fluctuations in the expanding phase 
(23). Current observations of large-scale struc- 
ture and fluctuations of the cosmic microwave 
background cannot distinguish between infla- 
tion and the cyclic model because both predict 
a nearly scale-invariant spectrum of adiabatic, 
gaussian density perturbations. 

Future measurements of gravitational waves 
may be able to distinguish between the two 
pictures (79). In inflation, where gravity is par- 
amount, quantum fluctuations in all light de- 
grees of freedom are subject to the same grav- 
itational effect described above. Hence, not 
only is there a nearly scale-invariant spectrum 
of energy density perturbations, but there is also 
a scale-invariant spectrum of gravitational 
waves. In the cyclic and ekpyrotic models, 
where the potential, rather than gravity, is the 
cause of the fluctuations, the only field that 
obtains a nearly scale-invariant spectrum is the 
one rolling down the potential, namely ¢, 
which only produces energy density fluctua- 
tions. The direct search for gravitational waves 
or the search for their indirect effect on the 
polarization of the cosmic microwave back- 


ground (29, 30) is the crucial test for distin- 
guishing inflation from the cyclic model. 


Cyclic Solution as Cosmic Attractor 


Not only do cyclic solutions exist for a range 
of potentials and parameters, but they are also 
attractors for a range of initial conditions. The 
cosmic acceleration caused by the positive- 
potential plateau plays the critical role here. 
For example, suppose the scalar field is jos- 
tled and stops at a slightly different maximal 
value on the plateau compared with the ex- 
actly cyclic solution. The same sequence of 
stages ensues. The scalar field is critically 
damped during the exponentially expanding 
phase. So by the time the field reaches stage 
3, where V = 0, it is rolling almost at the 
same rate as if it had started at d = 0, and 
memory of its initial position has been lost 
(11). The argument suggests that it is natural 
to expect dark energy and cosmic accelera- 
tion following matter domination in a cyclic 
universe, in accordance with what has been 
recently observed. 


Comparing Cyclic and Inflationary 
Models 


The cyclic and inflationary models have nu- 
merous conceptual differences in addition to 
those already described. Inflation requires 
two periods of cosmic acceleration: a hypo- 
thetical period of rapid expansion in the early 
universe and the observed current accelera- 
tion. The cyclic model requires only one pe- 
riod of acceleration per cycle. 

In the inflationary picture, most of the 
volume of the universe is completely unlike 
what we see. Even when inflation ends in one 
region, such as our own, it continues in oth- 
ers. Because of the superluminal expansion 
rate of the remaining inflating regions, they 
occupy most of the physical volume of the 
universe. Regions that have stopped inflating, 
such as our region of the universe, represent 
an infinitesimal fraction. By contrast, the cy- 
clic model is one in which the local universe 
is typical of the universe as a whole. All or 
almost all regions of the universe are under- 
going the same sequence of cosmic events, 
and most of the time is spent in the radiation-, 
matter-, and dark energy—dominated phases. 

In the production of perturbations, the 
inflationary mechanism relies on stretching 
modes whose wavelength is initially expo- 
nentially sub-Planckian to macroscopic 
scales. Quantum gravity effects in the initial 
state are highly uncertain, and inflationary 
predictions may therefore be highly sensitive 
to sub-Planckian physics. In contrast, pertur- 
bations in the cyclic model are generated 
when the modes have wavelengths of thou- 
sands of kilometers, using macroscopic phys- 
ics insensitive to quantum gravity effects. 

The cyclic model deals directly with the 
cosmic singularity, explaining it as a transition 
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from a contracting to an expanding phase. Al- 
though inflation does not address the cosmic 
singularity problem directly, it does rely implic- 
itly on the opposite assumption: that the big 
bang is the beginning of time and that the 
universe emerges in a rapidly expanding state. 
Inflating regions with high potential energy ex- 
pand more rapidly and dominate the universe. If 
there is a preexisting contracting phase, then the 
high—potential energy regions collapse and dis- 
appear before the expansion phase begins. 
String theory or, more generally, quantum grav- 
ity can play an important role in settling the 
nature of the singularity and, thereby, distin- 
guishing between the two assumptions. 

The cyclic model is a complete model of 
cosmic history, whereas inflation is only a 
theory of cosmic history following an as- 
sumed initial creation event. Hence, the cy- 
clic model has more explanatory and predic- 
tive power. For example, we have already 
emphasized how the cyclic model leads nat- 
urally to the prediction of quintessence and 
cosmic acceleration, explaining them as es- 
sential elements of an eternally repeating uni- 
verse. The cyclic model is also inflexible with 
regard to its prediction of no long-wavelength 
gravitational waves. 

Inflationary cosmology offers no informa- 
tion about the cosmological constant prob- 
lem. The cyclic model provides a fascinating 
new twist. Historically, the problem is as- 
sumed to mean that one must explain why the 
vacuum energy of the ground state is zero. In 
the cyclic model, the vacuum energy of the 
true ground state is not zero. It is negative and 
its magnitude is large, as is apparent from 
Fig. 1. However, we have shown that the 
universe does not reach the true ground state. 
Instead, it hovers above the ground state from 
cycle to cycle, bouncing from one side of the 
potential well to the other but spending most 
time on the positive-energy side. 

Reviewing the overall scenario and its im- 
plications, what is most remarkable is that the 
cyclic model can differ so much from the stan- 
dard picture in terms of the origin of space and 
time and the sequence of cosmic events that 
lead to our current universe. Yet, the model 
requires no more assumptions or tunings (and 
by some measures less) to match the current 
observations. It appears that we now have two 
disparate possibilities: a universe with a definite 
beginning and a universe that is made and 
remade forever. The ultimate arbiter will be 
nature. Measurements of gravitational waves 
and the properties of dark energy (//) can 
provide decisive ways to discriminate between 
the two pictures observationally. 
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Evidence of HIV-1 Adaptation 
to HLA-Restricted Immune 
Responses at a Population Level 
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Antigen-specific T cell immunity is HLA-restricted. Human immunodeficiency 
virus—type 1 (HIV-1) mutations that allow escape from host immune responses 
may therefore be HLA allele-specific. We analyzed HIV-1 reverse transcriptase 
sequences from a large HLA-diverse population of HIV-1—infected individuals. 
Polymorphisms in HIV-1 were most evident at sites of least functional or 
structural constraint and frequently were associated with particular host 
HLA class | alleles. Absence of polymorphism was also HLA allele-specific. 
At a population level, the degree of HLA-associated selection in viral se- 
quence was predictive of viral load. These results support a fundamental role 
for HLA-restricted immune responses in driving and shaping HIV-1 evolution 


in vivo. 


Selection of viral mutations associated with 
loss of antiviral cytotoxic T lymphocyte 
(CTL) responses has been described in hu- 
mans with acute and chronic HIV-1 infection 
(1), macaques infected with simian immuno- 
deficiency virus (SIV) (2, 3), and rhesus 
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monkeys challenged with simian-human im- 
munodeficiency virus (SHIV) after vaccina- 
tion (4). However, the full extent and im- 
portance of CTL escape mutation to HIV-1 
evolution remains to be established. CTL es- 
cape mutations occur at critical sites within 
HLA-restricted CTL epitopes where an 
amino acid substitution may abrogate 
epitope-HLA binding, reduce T cell receptor 
recognition, or generate antagonistic CTL re- 
sponses (/). Mutations that affect proteosome 
cleavage sites flanking CTL epitopes may 
also disrupt cellular processing of the epitope 
(5). The capacity of the virus to mutate at any 
amino acid residue is constrained, however, 
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